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ABSTRACT 

Gamma-ray bursts (GRBs) have been thought to originate from internal shocks 
that occur about 10^^ cm from a central site. The shells responsible for these shocks 
merge together and undergo an external shock at ~ 10^''^ cm, producing the afterglows. 
We include deceleration in our model of internal shocks and find that, for values of the 
Lorentz factor greater than 10^, deceleration is an effective catalyst for converting the 
bulk motion energy into radiation during the GRB phase. Substantial internal energy 
occurs because other shells run into the back of the first shell which has decelerated 
and because the first shell must energize the interstellar medium. Whereas internal 
shocks without deceleration are about 25% efficient, we can convert up to 85% of 
the bulk motion energy during the GRB phase. We demonstrate that the resulting 
time history can have three components. The first is due to internal shocks, but not 
those that involve the first shell. This component produces narrow peaks throughout 
the time history. The second is due to internal shocks involving the first shell, and 
it produces progressively wider and wider peaks but they tend to be hidden in a 
slowly varying background in the event. The third component is from energizing the 
interstellar medium. It is very smooth and may contribute mostly to a lower energy 
bandpass than the BATSE experiment. There have been claims of upper limits on the 
possible Lorentz factor because the deceleration must occur at greater radii than the 
internal shocks to avoid making progressively wider peaks. We do not find this to be 
the case, and the Lorentz factor can be much larger. 



Subject headings: gamma-rays: bursts 



1. INTRODUCTION 



Gamma-ray bursts (GRBs) are characterized by chaotic time histories which are often 
followed by x-ray, optical, and radio afterglows ( Metzger et al. 1997 , Costa et al. 199?! Frail et al, 
1997 ). The optical afterglows have shown redshifted absorption features which firmly established 
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a cosmological distance scale for the events. The distance implies that GRBs emit the order of 
10^^ to 10^^ erg (assuming isotropic emission). GRBs are also often characterized by emission up 



to 100 MeV with occasional reports of emission up to 10 Gev (Hurley et al. 1994). Given the 
photon density implied by the cosmological distances, photons above ~ 1 MeV would be destroyed 
by photon-photon attenuation if the emission is isotropic in our rest frame. Large relativistic 
bulk motion (Lorentz factors of ^ 100) allows for a much larger emitting surface combined with 
relativistic beaming that reduces the photon-photon attenuation ( [Fenimore, Epstein, &: Ho 1993 ). 



The high Lorentz factor plays a crucial role in virtually all models of GRBs. Originally, the 
prime suspect for the source of 10^^ erg was a neutron star-neutron star merger ( [Pacyznski 1986 ). 



However, such mergers were thought to occur on timescales (a few millisec) much shorter than 
GRB timescales (up to 10^ sec). [Meszaros fc Rees (1993) suggested that a relativistic shell would 



be formed by the initial release of energy. The shell could emit for a long time (10^ sec). If the 
shell is mostly moving directly at the observer, the shell stays close to the photons it emits such 
that they all arrive at a detector over a short period of time. If the shell moves with velocity v, 
then photons emitted over a duration At arrive at the detector compressed into a duration of only 
(c - v)At/c ^ At/{2T^) where T is the bulk Lorentz factor = (1 - ^2)"^/^ and [3 = v/c. The shell 
would emit due to the formation of an "external" shock when the shell decelerates by sweeping 
up the interstellar medium (ISM). In this explanation, density variations in the ISM cause the 
observed time structure. The deceleration is expected to occur at 

iJdec = b{pE^fl^Y^^''\m (1) 

where p is the ambient density (in cm~'^), is total energy (in erg) generated by the central 
site, and Lq is the initial Lorentz factor ( |Rees fc Meszaros 1994)) . For typical values such as p = 1 
cm~^, = 10^^ erg, and Lq = 100, the deceleration occurs at about 10^^ cm. The initial Lorentz 
factor is set by the baryon loading, that is 

£^0 = Tomoc^ (2) 

where mg is the mass of the shell, presumably carried by the baryons. 

An alternative explanation is that the central site produced a series of shells. Collisions 



between shells produce the gamma rays through internal shocks ( Rees fc Meszaros 1994 ). The 



faster shells catch up with the slower shells. The collision radius is roughly 

i?eoi = cT^AT (3) 



where AT is a typical time of variation in a GRB ( |Rees &: Meszaros 1994 ). For typical values 



such as AT = 0.1 to 1.0 s, R^oX is about 10^^ cm. The observed duration of the GRB is set by the 
duration of the activity at the central site. 



There are a series of arguments that indicate that the gamma-ray phase is not caused by 
external shocks. These arguments are all related to the fact that the size of the shell at the 
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deceleration radius is much larger than a causally connected region on the shell. For example, 
precursors and gaps require large causally disconnected regions to coordinate their activity 
( Fenimore, Madras, &: Nayakshin 1996| ). The observed variability implies that only a small fraction 
of the shell emits ( [Fenimore et al. 1996| , Sari & Piran 1997), and the average profile of many 
GRBs is inconsistent with that expected from a shell ( [Fenimore 1999| ). Finally, the constancy of 
the pulse width throughout the bursts indicates that we are not seeing a range of angles on a shell 



and that the shell is not decelerating ( Ramirez-Ruiz fc Fenimore 1999a ). A single decelerating 



shell would produce pulses that get progressively wider. 

However, external shocks have been very successful is explaining the x-ray, optical, and radio 
afterglows (see review by, e.g., piran (1999)[ ). These afterglows usually show power law decays and 
spectral variability expected from a decelerating shell. Thus, a general picture has formed where 
the central site produces multiple shells for tens of seconds. These shells collide, producing the 
gamma-ray phase by internal shocks and then merge into a single shell which interacts with the 
ISM to produce the afterglows ( [Sari &: Piran 1997| ). One-dimensional hydrodynamical calculations 
have reproduced many of the features of the spectral evolution ( Panaitescu fc Meszaros, 199^ ). 
There have been several detailed calculations of what is expected from the internal shocks. 



Mochkovitch, Maitia, k. Marques, (1995)[ , [Kobayashi, Piran, fc Sari (1997)[ , and paigne fc 



Mochkovitch (1998) used Monte Carlo calculations of internal shocks where the F, mass, time, and 
thickness of multiple shells are picked randomly to demonstrate that internal shocks can produce 
the variability seen in GRBs. By following their trajectories, it is determined when they collide. 
The duration of the resulting simulated GRBs is effectively the duration of the activity at the 
central site. The rise of each pulse depends on the time for a reverse shock to cross the shell and 



the fall depends on the curvature of the shell at the radius of interaction (Kobayashi, Piran, &; 



Sari (1997)). The radius of interaction depends on the Lorentz factors and the amount of time 
between the production of the shells at the central site. Indeed, the resulting time histories bear 
some similarities to the observed bursts. 

The strong optical emission discovery by the Robotic Optical Transient Search Experiment, 
ROTSE, ([Akerlof et al. 1999|) as served as an excellent test case for the external shock model. 



Sari &: Piran (1999) and Meszaros fc Rees (1999) fit forward and reverse external shocks and 



had excellent agreement with the time of the optical peak, the rise and fall times, the overall 
magnitude, and the break in the decay phase. 



The Role of F 



The Lorentz factor is not well determined observationally. The lack of apparent photon-photon 
attenuation up to ~ 100 MeV implies only a lower limit of ~ 100 ( Fenimore, Epstein, fc Ho 1993 ). 
The F determined by [Sari &: Piran (1999) for GRB990123 depended on some parameters adopted 
from Wijers &: Galama (1999) for GRB970508, but gave a similar low value of F = 200. However, 
there have been recent reports of possible TeV emission from GRBs ( [Leonor et al., 1999 ) implying 
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that r may be much larger. Also, if F is small, the efficiency of internal shocks is small, the order 
of 10%. 

The first shell starts to decelerate when it has swept up F"^ of its initial mass. Thus, larger 
F means the deceleration will occur at a smaller radius. But, a larger F means that the multiple 
shells will collide at a larger radius. Combining equations (|l]) and (^), internal shocks will occur at 
about the same radii as the deceleration if F^10~^*'(/?£'o)^''^. From the few observed redshifts, we 
know that GRBs have a rather large range of fluences, with typical values between 10^^ and 10^"^ 
erg for isotropic emission. We have little direct knowledge of the ambient density in the vicinity of 
a GRB, but it is reasonable to assume values of p between 0.1 and 10.0, so pE^ varies from 10^^ 
to 10^^ erg cm~^. For values of F of ~ 500 to 3000, the internal shocks will occur about the same 
place as the deceleration. 

Once the first shell decelerates (making an external shock), the rest of the shells will rapidly 



catch up to it resulting in rather efficient internal shocks. In previous models (e.g., Kobayashi et 



al. 1997| ), F was about 10^ to 10^ but deceleration was not included. It was assumed that the 



internal shocks would form at small radii and later the merged shell would suffer deceleration 
and an external shock. Perhaps a few straggler shells would catch up to the first shell after it 
decelerated and rejuvenate it during the afterglow phase ( [Panaitescu, Meszaros, &: Rees, 199^ ), 
but most of the gamma-rays were assumed to form at small radii relative to the external shock. 

We propose that the typical Lorentz factor is large enough such that the first shell decelerates 
before all of the multiple shells have a chance to collide. The deceleration is very rapid once it 
starts, effectively equivalent to slamming on the breaks. The rest of the shells catch up and collide 
with it. Since the efficiency of converting bulk motion to radiation in an internal shock depends 
on the difference of the colliding F's, the fact that the first shell is decelerating implies that the 
efficiency will be higher than in previous models. Thus, the collisions are internal shocks but the 
place and efficiency of the many of collision are caused by deceleration. 



3. Ingredients for a Model 



In the internal shock model, multiple shells are generated by an unspecified process at a 



central site. The parameters of our model will be similar to those of Mochkovitch, Maitia, & 
Marques, (1995), Kobayashi et al. (1997) and Daigne &: Mochkovitch (1998), including the time 
the i-th shell was generated (to«)) the initial width of the shell (/j), the minimum and maximum 
initial Lorentz factor (Fmin, Fmax)) and the initial energy {E,j). Kobayashi et al. (1997) allowed for 
selecting the initial mass, energy, or density. However, all three gave similar results and we will 
restrict ourselves to selecting the energy. The initial rrii is then found from equation (§). Since 



Kobayashi et al. (1997) presented unitless intensities, it was unnecessary for them to specify Ei. 



The bulk energy is necessary to set the deceleration. The peak energy can be estimated from 
bursts with observed redshifts. GRB970508 had a peak luminosity, L, of ~ 3 x 10^^ erg s~^. Other 
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GRBs have shown extreme redshifts (e.g., |Kulkarni et al. 19991) , implying L ~ 2 x 10^^ erg s~^. 
We will uniformly select Ei between -Emin (=10'^^ erg s~^) and -Emaxi and vary -Emax from 10^^ to 

IQ53.5 gj,g 

We randomly select tot+i — toi from a Poisson distribution based on the rate of peak 
occurrence. Thus, we specific the duration of the activity at the central site (Tjur) and the 
expected number of peaks (N) such that the actual number of peaks is random. Since they were 



presenting results in unitless time, Kobayashi et al. (1997) set the burst duration, shell separation 



and the number of peaks to be constants. These difi^erences are not important when there is no 
deceleration. We will use parameters that are roughly equivalent to [Kobayashi et al. (1997) , that 



is, = 85, Tdur = 60 s, and k = 0.2 s. Burst often have gaps which implies that the activity at 
the central site can turn off for a while. To demonstrate the effects of turning off the central site, 
we impose a gap in the activity between T = 20 and T = 33 s. 

Until they collide with the first shell or each other, the motion of the every shell except the 
first is constant, Ri{t) = cf5i{t — toi). If Tj is greater than Tj, the two shells will collide at time tij 
when Riitij) = Rj{tij), which we call the collision radius {=Rc), and it occurs at 

« 3 

where Afoij = ^oi ~ ^Oj- The resulting pulse arrives at a detector at the relative time of arrive 

?toa — ^ij ~ Rcl C — toi + =^ Atgij . (5) 

« ~ j 

Thus, the relative time of arrival at a detector will have a close one-to-one relationship with the 
time the shell was created (i.e., t^i). 

In order to conserve both momentum and energy when shells collide, some of the bulk energy 
must be converted to internal energy which will be radiated away. If -Erad is the generated internal 
energy, conservation of energy dictates that 

miTj + mjFj = [rriij + ^^'^ij (6) 

where Fjj is the Lorentz factor of the resulting shell and the resulting mass is rriij = rrii + mj. 
Conservation of momentum gives: 

rriiPiTi + rrijPjTj = [niij + -^IftjFjj (7) 

where, as usual, the F terms are related to the /? terms as F = (1 — P'^)^^^'^. The post collision /3 is 

ruiTi + irijTj 
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which has the approximate solution (Kobayashi et al. 1997) 

iriiTi + m,T 



(9) 



The first shell is decelerated by the ISM. We use equations ^ and ^ with Tj = 1 and mj equal 
to the mass swept up during the time step to determine the velocity of the first shell as a function 
of time. 

The colliding shells make a peak in the gamma-ray time history at relative time tij — Rc/c. 
When two shells collide, forward and reverse shocks traverse the shells. If the internal energy is 
promptly converted into radiation, the merged shell emits for about the time that it takes for the 



[1 (see 


Kobayashi et al. 1997" 


|Sari & 


Piran (1995) 





and 



1 + 2Tij/Ti 

2 + Tij/Ti 

1 + 2Tij/rj 

2 + T,,/Tj 



1/2 



1/2 



(10) 



(11) 



4. Pulse Shape 

A shell that coasts without emitting photons and then emits for a short period of time 
produces a pulse with a rise time related to the time the shell emits and a decay dominated by 
curvature effects ( Fenimore et al. (1996)| ). In the internal shock model, the shell emits for the time 



it takes the reverse shock to cross the shell that is catching up, that is ( [Kobayashi et al. 1997| ), 

At,,oss = lj/{Pj - Prs) ■ (12) 

The time of arrival at a detector (relative to the start of the pulse) of photons generated at 
angle from the line of sight is 

T{e) = Rc{l-cose)/c (13) 

(Note in our previous papers, it was more convenient to measure time from when the shell left the 
central site; this is not used here because the shell does not move at a constant speed.) At angle 
9, the Doppler factor. A, is Fjj(l — Pij cos9). At time T in the pulse, the A factor is 

Rc + 2F2. cT 

A(T) = \ ^ (14) 

To calculate the observed pulse shape, one needs to combine the Doppler beaming with the volume 
of material that can contribute at time T. Following the method in Rummer fc Fenimore (1998) 
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the resulting pulse shape is 

v{T) = ifr<o 

= i'- ; — ^ ^-TT^ — if < 2rlT < Atcross (15) 

where -0 is a constant and T is measured from the start of the pulse. 

The cooling is very rapid so the internal energy generated by the collision, -Erad) is immediately 
turned into photons. An observer, using a detector such as BATSE, sees the fraction that is in the 
BATSE bandpass of 50 to 300 KeV, /batse- Since we do not understand exactly how the internal 
energy is distributed, we cannot predict /batse- However, GRB often have a "Band" spectral 
shape with a = — 1, /? = —2.5, and i^pcak = 250 KeV ( [Band et al. 1993 ). If that shape is valid over 



the entire range of emission, /batse is ~ 0.37. We generate simulated time histories as the sum of 
pulses with the shape from equation (|l5|) and integrated fluence of /BATSE-E-rad- We generate the 
time history with 0.064 s samples (to mimic BATSE) and then find the peak emitted luminosity 
in 0.256 s (= -^256)- Ignoring cosmological redshift effects, the BATSE catalog value of P256 should 
be related to ^256- 



5. A Typical Simulation 

To summarize our model, we have eight parameters: the duration of the activity (Tdur), the 
rate of explosions at the central site (A^/Tdur), the range of Lorentz factors (rmin,rmax)) the range 
of energy release at the central site (-Emim -E'max), the ISM density (/)), and the range of initial 
thicknesses (0 to /). To be comparable to [Kobayashi et al. (1997)| , we will set A'" = 85, T^uv = 60 
s, and I = 0.2 s. [Kobayashi et al. (1997) parameterized much of their results based on rmax/rmin 



since the overall efficiency of the conversion of bulk energy to radiation was primarily dependent 
on that parameter. With deceleration, we have found that the efficiency depends mostly on Fmax, 
so we have set Tmin to the minimum required for the high energy emission (100), and varied Fmax 
from 10^'^ to lO'^'^. Kobayashi et al. (1997)[ had no analog to -Emin, -E'max) and p. Since we selected 



E uniformly between ii^min and -Emaxi -E-min is not important as long as it is much less than -Emax- 
We set ii^min to 10^^ erg. For p we have used 1 cm~^. 

Figure ^ is a typical simulation -Emax = 3.6 x 10^^ erg, F^ax = 3.2 x 10^. In Figure |l]a, p is 
zero, so there is no deceleration. About 1.5 x 10^^ erg (assuming isotropy) were released at the 
central site in 81 shells. The burst duration at the observer is approximately the duration of the 
activity at the central site. About 25% of the bulk energy was received by the observer in the 
period Tdur- Figure Qb is the same simulation (i.e., same set of random numbers), but includes 
deceleration of the first shell in an ISM with p = 1 cm~^. Both simulations appear similar because 
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both reflect the activity of the central engine (see eq. ||5|). The dotted line is the contribution to 
the time history from collisions with the first shell. It tends to add a DC level with a few wide 
peaks but it raises the fraction of the bulk energy converted to radiation to 45%. Some of the 
radiation will arrive after T^ur because curvature will delay it. 

Figure ^ gives the Lorentz factor for the first shell in Figure |l|b. Given the high value of 
Fmaxj it quickly decelerates but other shells collide with it, giving it a boost and maintaining a 
large F for most of the burst. The deceleration occurs because the first shell collides with the ISM. 
The resulting internal energy must also radiate away. In Figure ^3 we show the contribution to 
the time history of Figure |l| from the internal energy from the deceleration if it radiates in the 
BATSE bandpass. It tends to be smooth and would fill in gaps if it had a /batse similar to that 
from the internal shocks. We define -Edecjciur to be the internal energy generated by the collision 
of the first shell with the ISM that would arrive at the detector with Ttoa < ^dur- For the case in 
Figure |l|, -Edec,dur is 38% of the bulk motion energy. 



6. Efficiency of Converting Bulk Energy 

The efficiency of the conversion is an important constraint. Although the time histories imply 
that GRBs are central engines with internal shocks, internal shocks usually do not convert most of 
the bulk motion into energy (e.g., ^25%, [Kobayashi et al. 1997] ). Observationally, the afterglows 
only account for a small percentage of the energy so it is not clear where most of the energy goes. 
The efficiency for an individual collision can be found from the initial and final bulk energies: 

e = 1 ^ (16) 

mjFj + mjVj 

If there is no deceleration, the shells will collide until the remaining shells are ordered with 
decreasing value of the Lorentz factors. Let n be the number of remaining shells. The overall 
efficiency depends on how much energy remains in un-collided shells: 



_ ij=0 ^ij^ i] 

l^i=o rriiL oi 



When deceleration occurs, n is 1, another reason why our model will give higher efficiency than 
previous models. 

To study the effects of deceleration we have generated sets of 128 bursts, under a variety 
of conditions. Figure H shows the average efficiency as a function of Fmax- The curves labeled 



"No Deceleration" is effectively the same result as Kobayashi et al. (1997) . The curves labeled 
"Deceleration, IS" includes an ISM with p = 1 cm^'^. We ran models for a range of E'max 
(maximum energy per shell) and interpolated the results to find the efficiency at three values of 
the peak L256 in the time histories: 3 x lO^'', 3 x 10^^, and 3 x 10^^ erg s^^ (the solid, dotted, and 
dashed lines, respectively). Figure ^ shows the corresponding average radii for the internal shocks 
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that produces pulses that arrive with Ttoa < T^ur, that is, during the GRB phase. (Including all 
internal shocks would produce a misleading result when there is little deceleration because a few 
stragglers would finally collide at radii orders of magnitude larger.) The curves labeled "ES" are 
the average radii at which the Lorentz factor of the first shell is reduced by half. Once the first 
shell starts to decelerate, a fair number of the shells collide with it, raising the average amount of 
the bulk energy which is released in internal shocks. These collisions are more efficient since there 
is a greater disparity between the T factors. For large vales of Tmax; the efficiency rises to 40%. 

The curves labeled "Deceleration, IS" in Figure |3| are based on the ratio of the internal energy 
generated by collisions between shells (including the first shell) to the total generated bulk motion 
energy. It does not include the bulk motion energy lost to energizing the ISM. Eventually, all of 
the bulk motion energy is lost to the ISM. In previous models it was assumed that this was far 
from where the internal shocks occur. In the curves labeled "Deceleration, IS+ES", we include 
-E'dec,rad in the efficiency, that is, the internal energy from collision with the ISM whose photons 
would start to arrive during the burst. For large values of Fmax) nearly 85% of the bulk motion 
energy is lost during the GRB phase. 



7. DISCUSSION 

Internal shocks are capable of producing the variability that is the signature of GRBs 
(Kobayashi et al. 1997). However, it has been believed that internal shocks are inefficient. 



converting only ^25% of the bulk motion energy into radiation. Since the afterglows only account 
for a few percent of the radiated energy, it has been unclear where most of the energy goes. 

In this paper, for the first time, deceleration of the first shell is included in an internal shock 
model. For Fmax^lO^, there are two ways that deceleration is an important catalyst for converting 
bulk motion into radiation. First, the deceleration occurs because the bulk motion must energize 
the ISM that it runs into. Much of the energy to energize the ISM goes into internal energy. This 
is rather effective because the bulk motion energy of the first shell is FjMjC^ where the mass grows 
as other shells run into the first shell. If deceleration causes the F of the first shell to drop by 50%, 
nearly 50% of the bulk motion energy will be used. Second, the rapid deceleration causes shells 
to plow into the back of the first shell. The efficiency for converting bulk motion (for equal mass 
shells) scales as 1 — (Fj/Fj)^/^ if Fj is much larger than Tj, as is the case when the j-th shell is 
decelerating. These two effects combine to release up to 85% of the bulk motion energy while T^a 
is less than T^un that is, during the GRB phase. Although the bulk motion might be effectively 
converted to internal energy, the resulting electron distribution is difficult to predict, so we do not 
know how much of it will occur in the BATSE bandpass (i.e., /batse is uncertain). 

Thus, one can identify three types of contributions to the time history, each with a different 
character. The internal shocks that do not involve the first shell, internal shocks involving the first 
shell, and the external shock produced as the first shell decelerates. 
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The internal shocks that do not involve the first shell are characterized by narrow pulses, and 
have nearly constant width throughout the time history (see, for example. Fig. |^a). The Ttoa 
for these pulses is dominated by the time the shells were produced at the central site (following 
eq.[|[). They occur at a similar radius from the central site. If Fj is selected randomly between a 
small Fniin and Fmax, many pulses form with a similar Lorentz factor: F^ ~ (FjFj)^/^ ~ Fniax/2. 
The pulse shape depends mostly on Rc and Fjj (eq. |0]), so they are quite similar. 

The internal shocks involving the first shell occur at ever increasing radii with a generally 
decreasing Lorentz factor. Thus, equation (|l^) produces peaks that are wider and wider (see the 
dotted curve in Fig. In our previous papers, we argued that the time history could not arise 
from a single shell because the pulses did not get wider and wider. This argument is still valid 
and this paper shows how multiple shells can produce narrow peaks throughout the event in the 
presence of wider and wider pulses from a single shell. Indeed, a recent analysis of 387 pulses 
in 28 BATSE GRBs shows that the most intense pulses in a burst have nearly identical widths 
throughout the burst, but the weak pulses show a trend to become wider as the burst progresses 
(Ramirez- Ruiz &: Fenimore 1999b). This is precisely what is seen in simulations when deceleration 



is included. 

There have been claims of upper limits on the possible Lorentz factor because the deceleration 



must occur at greater radii than the internal shocks ( [Lazzati, Ghisellini, fc Celotti 1999]) to avoid 
making progressively wider peaks. We do not find this to be the case, and the Lorentz factor can 
be much larger, allowing more of the bulk motion energy to be released during the GRB phase. 

The third type of contribution arises from the external shock as the first shell energizes the 
ISM. This is a smooth component with some variation as the Lorentz factor increases due to 
collisions with faster shells and decreases due to deceleration (see Fig. ^3). Previous external 
shock models (e.g., Dermer &: Mitman 1999| ) have assumed the shock interacts with ISM clouds 



that are much smaller that the size of the shell. This was necessary to produce the temporal 
variability. We do not assume any structure in the ISM so the contribution from this component 
is quite smooth. 

Figure |5| shows BATSE time histories that have the characteristics of our simulations. Figure 
^a is BATSE burst 2831. It has many narrow peaks throughout the time history but also gaps that 
go back to background. The presence of gaps implies little deceleration of the first shell because 
the gaps are not filled in. The gaps would occur because the central site turns on and off. We note 
that this burst is the record holder for the highest energy photons, 18 GeV ( Hurley et al. 199^ . 



Figure ^b is BATSE burst 2329, and initially, it shows narrow peaks, but then broader peaks, 
on top of a slower varying level. The structure on the rise in burst 2329 is statistically significant. 
The slower varying function and widening pulses implies substantial deceleration of the first shell. 
In other bursts (e.g., BATSE burst 130), there are gaps where the last pulse before the gap has 
a slow decay, similar to that seen in Figure |T|b. However, these gaps can go all the way down to 
the background. Such bursts seem to show the signature of internal shocks on a decelerating first 
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shell but not the contribution from the energization of the ISM. Since we do not understand the 
mechanism by which the internal energy is distributed, the /batse values associated with each 
component might be different and the emission might appear in different bandpasses. Perhaps the 
/batse for energizing the ISM is small and its internal energy is radiated at lower energy such as 
the x-ray excesses reported by Preece et al. (1995). 



Apparently, some burst involve deceleration and some do not. The Lorentz factor required 
to have deceleration depends very weakly on p and £"0 (i-e., {pE^Y^'^). Thus it seems more likely 
that intrinsic variations in Fmax might be the reason why some bursts show more deceleration 
than others. 

If the prompt emission is caused by the first shell, as suggested by the analysis of GRB990123 



(Sari &: Piran 1997), we would expect events with a slowly varying component to be more likely 



to slow prompt, bright optical emission or early afterglows. 

In summary, it is possible to convert a large fraction (~ 85%) of the bulk motion energy into 
radiation during the gamma-ray burst phase with internal shocks if deceleration of the first shell 
is account for and the Lorentz factor is ^10^. 
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Fig. 1. — Simulated GRB time histories, (a) A simulated GRB time history from internal shocks 
with no deceleration. The Lorentz factors were chosen uniformly between 10^ and 10^'^ and the 
peak luminosity (found over 256 ms) is ^256 = 2 x 10^^ erg s~^. Shells were randomly created at 
the central site for 60 sec except for a 13 s period, thus producing a gap in the simulated GRB 
time history. About 25% of the bulk motion was converted into radiation by the internal shocks, 
(b) Simulated GRB time history from internal shocks including deceleration of the first shell due 
to an ISM with 1 particle cm~'^. The same random numbers as in (a) were used. The lower curve 
is the contribution to the time history from collisions with the decelerating shell. Note that the 
there are still many narrow spikes throughout the event. There is a slower decay during the gap. 
About 45% of the bulk motion was converted into radiation by the internal shocks. 

Fig. 2. — (a) The Lorentz factor of the decelerating shell in Figure |l| as a function of when photons 
would arrive from it. The Lorentz factor is fairly constant for much of the burst because shells 
continue to collide with it. (b) The expected time history from the bulk motion energy that is 
converted into internal energy by the collision of the first shell with the ISM. It is smooth and 
would fill in gaps in the time history. 

Fig. 3. — The efficiency of converting bulk motion into internal energy as a function of the maximum 
Lorentz factor. If there is no deceleration (ISM density = 0), the maximum efficiency is about 25%. 
If deceleration is included, the internal shocks convert up to 45% of the bulk motion to internal 
energy. The deceleration is caused by an external shock that sweeps up and energizes the ISM. The 
curve labeled "Deceleration, IS+ES" gives the fraction of the total original bulk energy that is lost 
from the shells (some into internal shocks, some to energized the ISM). Up to 85% of the original 
energy is used during the GRB phase if the Lorentz factor is as large as 3 x 10^. 

Fig. 4. — The radii for deceleration and internal shocks. The solid, dotted, and dashed lines are for 
increasing values of L256- The curves labeled "ES" are the average radii when the Lorentz factor 
of the first shell first drops to half its original value. The curves labeled "IS" are the average radii 
for internal shocks that produce peaks during the GRB phase. For Lorentz factors greater than 
^ 3000, many of the internal shocks are caused by shells running into the back of the first shell 
which is undergoing deceleration. 

Fig. 5. — The time history of long BATSE GRBs. (a)BATSE burst 2831 which is consistent with 
internal shocks with very little deceleration. Note that the gaps are near background and the decay 
into the gaps is sharp. (b)BATSE burst 2329 which is consistent with substantial deceleration. 
Note that the peaks tend to get progressively wider and there is an underlying smooth component. 
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